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1. Introduction 


In several recent studies, asymptotic methods have been used success- 
fully for the derivation of rat.onal approximations which describe the inter- 
action of a turbulent boundary layer and a weak, stationary, normal shock 
wave. It appears that correct limiting forms of the equations can be deter- 
mined, that numerical or analytical solutions to these equations are ob- 
tained easily enough to be of practical interest, and that numerical accu- 
racy may be adequate for important parameter ranges. In the limiting 
case to be considered here, still for an unseparated boundary layer, 
the shock wave extends close to the wall, the upstream influence is 
small, and analytical solutions can be obtained for most of the flow 
field. Pressure distributions are derived in Part I; the wall shear 
stress and the possibility of predicting separation will be discussed 
in Part II. 

In many transonic flows of interest, there occurs a shock wave which, 
in an inviscid-flow approximation, is norrhal to a solid boundary, at values 
of the Reynolds number large enough that the boundary layer along the wall 
is fully turbulent. Since the strength of the shock wave must decrease to 
zero in the supersonic part of the boundary layer, there can be no discon- 
tinuity in the pressure at the wall. It is observed that the shock wave 
becomes slightly curved and is displaced slightly in the upstream direc- 
tion. As the Mach number upstream is increased, still below the value 
required for separation, the shock wave extends further into the boun- 
dary layer; experimental results [1,2] show an initially rapid rise 
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in the wall pressure, followed by a gradual decrease in the pressure 
gradient over a distance several times larger than the boundary layer 
thickness. 

Asymptotic descriptions of these flows, in the limit of infinite Reynolds 
number, have been discussed in References [3] through [lO]; in particular. 
Ref. [8j contains the first steps of the present work. In each of these 
studies, the representation of the undisturbed boundary layer in terms of a 
velocity-defect layer and a wall layer [ll, 12, 13] is regarded as providing an 
asymptotic description as the Reynolds number tends to infinity [l4- 18]. 

The pressure gradient in the boundary layer is large near the shock wave, 
and consequently the forces resulting from changes in the Reynolds stresses 
are of higher order than terms retained, in most of the boundary layer. 

Thus, as for laminar flow [19, 20, 2l] , an asymptotic description of the 
changes in the mean flow can be obtained with the use of inviscid-flow equa- 
tions for most of the boundary layer. 

The form of the velocity profile, however, implies two important 
differences from the laminar-flow case. First, for an unseparated turbu- 
lent boundary layer the wall layer is extremely thin, and the displacement 
effect resulting from deceleration of fluid close to the wall remains very 
small, even in a large pressure gradient. Thus, if the undisturbed velocity 
profile is known outside the wall layer, an approximation to the pressure 
can be found without knowledge of the flow details near the wall and there- 
fore without any further assumption about the nature of the turbulent 
stresses. Second, for a slightly supersonic external flow the sonic line is 
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located at axi arbitrary position (outside the wall layer) in the undisturbed 
boundary layer, depending on the relative sizes of the nondlmensional fric- 
tion velocity and the nondimensional difference between the fluid velocity 
and the critical sound speed in the external flow. As the Reynolds number 
tends to infinity, one can then study three cases, such that the ratio of 
these parameters tends to infinity, remains constant, or approaches zero. 

Adamson and Feo [3] considered an incident oblique shock wave in a 
flow with velocity only slightly greater than the sound speed, such that the 
sonic line is located very close to the edge of the boundary layer. The 
corresponding asymptotic formulation was shown to lead to a local- inter- 
action problem requiring solution of the transonic small-disturbance equa- 
tions for the local perturbations in the external flow, expressed in appro- 
priately scaled variables. The influence of the boundary layer is repre- 
sented on this scale through an effective wall boundary condition specifying 
a linear relationship between the streamline slope and the pressure gradi- 
ent. Melnik and Grossman [4] studied a normal shock wave having strength, 
as measured by the nondimensional pressure jump, of the same order as 
the friction velocity, so that in the limit the sonic line is at an arbitrary 
location in the boundary layer. Numerical solutions of the transonic small- 
disturbance equations were obtained for perturbations in the defect portion 
of the boundary layer and in the neighboring external flow. Changes in the 
wall layer were also discussed in each of these papers. Melnik and Grossman 
later [5, 6] obtained additional numerical solutions for axisymmetric pipe 
flow. At higher upstream speeds, which might be characterized as "high 
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transonic speeds," the shock wave is stronger bat the boundary layer can 
remain unseparated. For this case, a first approximation for the flow per- 
turbations outside the wall layer was given by Adamson and Messiter [8] . 

The shock-wave strength, although still small, was taken to be large in com- 
parison with the nondlmensional friction velocity, so that in the undisturbed 
boundary layer the distance from the sonic line to the wall is much smaller 
than the boundary-layer thickness. The corresponding problem has also 
been discussed for an incident oblique shock wave [7, 9] • A brief prelimi* 
nary description of some of the present results was given in Ref. [lO] ; a 
few details have since been modified. 

In the present work, analytical solutions are derived which incorporate 
additional physical effects as higher-order terms for the cas^ first dis- 
cussed ratiier briefly in Ref. [8] , when the sonic line is very close to the 
wall. The functional form used for the undisturbed velocity profile is des- 
cribed in Section 2, to indicate how various parameters will be calculated 
for later comparison with experiment. The basic solutions for the pressure 
distribution are derived in Section 3. In Section 4 corrections are added 
for flow along a wall having longitudinal curvature and for flow in a circular 
pipe, and comparisons with available experimental data are shown. 
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2. Undistarbed Velocity Profile 


Nondimensional rectangular coordinates X and Y are measured along 
and normal to the wall, respectively, with Y = 0 at the wall and X ~ 0 at 
some point on the shock wave, e« g. , at the intersection of the shock wave 
with the edge of the boundary layer as defined below. The reference length 
is a geometric length such as the length of the boundary layer from a lead- 
ing edge up to the shock wave. The nondimensional mean-velocity compo- 
nents U and V, referred te the critical sound speed in the external flow, 
are in the X and Y direcUcns respectively, cind the term p' V'/p has been 
included in V. Here primei' denote fluctuations about the mean, and p ' V' 
denotes an average value. The nondimensional mean pressure P, density 
p, temperature T, and viscosity coefficient p. are referred to the criti- 
cal values of pressure, density, and temperature, and the corresponding 
viscosity coefficient, in the flow just outside the boundary layer and ahead of 
the shock wave. The sum of the nondimensional Reynolds stress and viscous 
stress,' in the boundary-layer approximation, is denoted by t, and has been 
made nondimensional with twice the dynamic pressure, in terms of the same 
reference quantities. For later convenience the friction velocity u^ is made 
nondimensional using the external-flow density; 
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U =1+6 
e 


( 2 . 1 ) 


where the subscripts e and w indicate values in the external flow and at the 
wall, respectively, and c^ is the undisturbed value of the skin friction 
coefficient, referred as usual to the dynamic pressure in the external flow. 
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The nondimens ional difference between the fluid velocity and the critical 
sound speed in the external flow is c , and in the present case u^ « e « 1. 
For simplicity, an adiabatic wall is assumed and the total enthalpy is taken 
to be uniform. The ratio of specific heats is \ is constant. 


As in references cited above, it is assumed that the undisturbed 
boundary layer can be described asymptotically in terms of a velocity-de- 
fect layer and a wall layer. The defect layer occupies most of the boundary 
layer, and its thickness is taken equal to a boundary-layer thickness 6. 

The velocity differs from the external-flow velocity by an amount of order 

2 

u , the shear stress is t = 0(u ), and the layer thickness is 6 = 0(u ). 

T T T 

The much thinner wall layer has thickness denoted by 6, and the velocity 
there is small, of order u^. Coordinates measured in terms of these non- 
dimensional thicknesses are defined by 


y = 


Y 
6 * 


~ Y 

y = r 


6 = O(u^) 

~ W.2 e 1 

•^e e T 


( 2 . 2 ) 


(2.3) 


where 6 « 6 , and 6 has been set equal to the ratio of the nondimens ional 
local kinematic viscosity and a friction velocity u (T / T )^^^ = (t /p 
based on the density at the wall. The Reynolds number Re is based on the 
geometric reference length and the undisturbed external-flow velocity and 
kinematic viscosity; all parameters are understood to be evaluated imme- 
diately upstream of the shock wave. 
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The velocity in the undisturbed boundary layer just ahead of the 
shock wave is expressed in the defect layer in terms of y and in the wall 
layer in terms of y, as follows: 


"a ~ “e 


y = 0(1) 


U.. ~ ./ TJ*'' * y = 0(1) 


u T w e 


(2.4) 
(2. 5) 


The form of the profile is shown in Fig. 1 for u^ « e « 1. Equations 
(2.4) and (2, 5) are [l3] , respectively, the ’’law of the wake" and the "law 
of the wall," written here for a compressible boundary layer, and are 

/Sit 

taken to be asymptotic representations valid as u^ -► 0, with y and y held 

fixed respectively. Throughout most of the analysis also e 0 such that 

u^/ € 0. In the wall layer the Reynolds stress and the viscous stress 

2 

are both of the same order as the wall shear stress t = 0(u ). Since 

w T 

/St 

Y = 0(6) is extremely small, the momentum equation gives t ~ As 

/St » tst 

y = y/ 6 00 , the viscous stress becomes extremely small, while x re- 

mains equal to in the limit, provided that also y = Y/ 6 -► 0. The mix- 

2 t ^2 

ing length approximation K p (ydU^/ dy) = x^ + . . . is introduced here for 
y « 1 and y X, where k is the von Kftrman constant, taken equal to 

0, 41. For a perfect gas with uniform total enthalpy, p T = and 

112 ~ 

T = + 1) - - 1)U , Integration gives, for y « 1 and y » 1, 


1 , 

.2 


= rsin{r-‘(T^Ty u^(/c'*toy + c)} 


( 2 . 6 ) 


l/2 1/2 

where c = constant and F = (y + 1) ' /(y - 1) ' . This is van Driest’ s 
[22] result, with the added simplifying assumption of uniform total enthalpy. 
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Expansions of Eqn, (2. 6) for -*■ 1 + c and for -*■ 0 should agree, 
respectively, with expansions of the defect-layer velocity (2. 4) as y -* 0 and 
of the wall'layer velocity (2. 5) as y oo. For -► 1 + g and -► 0, 
respectively, Eqn. (2.6) gives 

~ 1 + € + (u ^K)(ln y “ 2II ) (2.7) 

^u ~ In y + c) (2. 8) 

where II is Coles* [13] profile parameter; c « 5.0 and, for zero 

pressure gradient, II =s 0. 5 or perhaps a little larger. Since y = (6/6)y, 
comparison of £qns« (2. 6) and (2. 7) gives 

X ■ Vfi A) * ) - U^(2nx + c) (2,9) 

where U.(c) = F sin ^(F ). The expansions (2.7) and (2.8) require, 

X G 

respectively, y -► 0 slowly and y -► oo slowly as u^ -*■ 0; since 

= 0(1 /in Re), from Eqn. (2.9)» one might take, e.g. , y = 0(u^ 

and y = 0(u^”) as u^ 0, where m > 0 and n > 0. A difference from the 

incompressible case arises because Eqn. (2.8) with y = (6/6)y does not 

^ 

agree with Eqn. (2*7). That is, the expansion as y -* oo of the wall^layer 
solution does not agree with the expansion as y -► 0 of the defect-layer solu- 
tion. Thus these solutions have no common domain of validity and cannot 
be matched. This type of problem has been discussed in detail by Lager - 
Strom eind Casten [23] , with a model example related to flow at low Reynolds 
number. In the present case, the density has different values for y = 0(1) 
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and for y = 0(1), and the difficulty is resolved by use' of the solution (2. 6) 

» 

for 6 « Y « 6 ; this feature was also noted by Adamson and Feo [3] and 
by Melnik and Grossman [4] . 

The defect layer, where y = 0(1), has nearly constant density and is 
described in a first approximation by incompressible -flow equations. The 
domain of validity of Eqn. (2. 6) can be made to include y - 0(1) if 
K ^(In y - 2n) is replaced by ^^^(y), where is the same 

■ .f ■ . . 

for incompressible flow. Then 

= rsin{sta'V'\) + 

Expansion for u^ -*■ 0 gives Eqn. (2,4) if y is held fixed, Eqn. (2,7) if 
y -► 0 sufficiently slowly that also u_^ In y -*• 0, and Eqn. (2.6) if 
y = (6 /6)y -► 0 more rapidly, such that u^ In y is held fixed. The use of 
Eqn. (2. 10) was suggested by Maise and McDonald [24] , who showed that 
this assumed profile permits good correlation with experimental data for 
adiabatic flat-plate boundary layers. Their interpretation of Eqn. (2. 10) 
notes that a transformed velocity F sin ^(F ^U^) is predicted to have the 
incompressible form U. (c ) + (T /T ) ' u o-^(y) everywhere outside the 

X © T vX 

wall layer. 

A second relation between 6 and u_ for 9 P/8 X » 0 can be found with 

T 

the help of the von Karman integral of the momentum equation, following a 
derivation similar to that for incompressible flow given, e. g. , by 
Cebeci and Smith [25] . The result is, to second order in u 
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m. 

T 


-6 = {T u. ^ ^ “oi'^y^ ^2 


The positive constant nl^ is defined by 


oo 


m^ = •* ® «oi^y)dy 


( 2 . 12 ) 


and occurs in aT^other context in the following section. For analytical 
purposes, the function iiQj^(y) is represented in Coles' [l3] form' 


Uoj^(y) = K ^ In y ^(1 + cos Try) (2.13) 

for 0 < y < 1, with ® for y > 1. 

In the derivation which follows, the boundary-layer thickness is taken 
as one of two Important characteristic lengths. The other leng& is the dis- 
tance from the wall to the sonic line in the undisturbed boundary layer, de- 
noted in nondimen sional form by 6^. Substituting Eqn, (2. 13) in Eqn, (2. 10), 
setting y = 6jj,/6, and expanding for 6^ /6 ^ 0 gives 


u^x '^ln(6/6 ^) ^ ) - U.(0)] - 2Jl ^ ‘^u^ (2. 14) 

As e 0, In (6/6^) ~ k e Jl - (y - l)e/4 + . . .] - 2ll ; thus 6jj./6 0 

if u^/e -»• 0, An alternate form of the velocity profile (2. 6) in terms of a 
coordinate y = 

^ r 8 in{sin“^(r’^ ^ y*> <2. 15) 

rM • 

for y *♦ 0 and y -*• oo. 
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3. Interaction Along a Plane Wall 


As “*• 0, the orders of magnitude of the mecin pressure gradient and 
fluid acceleration near the shock wave are larger than In the undisturbed 
boundary layer. The Reynoldsrstress transport equations can be used to 
show that in most of the boundary layer the contributions to the mean forces 
resulting from changes in the turbulent stresses are sufficiently small» in 
comparison with the pressure and inertia terms, that they may be neglected 
as u^ 0, not only in a first approximation but also in the calculation of 
some higher-order terms. Correct asymptotic representations of the mean 
velocity and pressure perturbations can therefore be derived using inviscid- 
flow equations. Also, as noted at the end of this section, displacement 
effects resulting from flow changes very close to the wall are ext'remely 
small, and so the largest terms in the solution for V should approach 
zero as the distance from the wall decreases. 

In the equations which follow, all laminar and turbulent stresses are 
neglected, as are the entropy changes across the shockwave; order-of- 
magnitude estimates given at the end of this section show that the neglected 
terms are in fact of higher order than any of the terms retained. The equa- 
tions describing the fluid motion can then be written in the following form; 

a^ div q = q • V ^ (3.1) 

a^ = |- (\ + 1) - - l)q^ (3. 2) 

pj • V q = -y V P (3.3) 


U 



-* 1/2 

Hereq, q, and a = (P/p) ' are, respectively, the velocity vector, the 
magnitude of the velocity, and the sound speed, all nondimensional with the 
critical sound speed in the external flow jurt ahead of the shockwave. Tlie 
gradient and divergence operators imply differentiation with respect to the 
nondimensional variables X and Y. Crocco ' s theorem, simplified by tiie 
assumption of uniform total enthalpy, is 


n X q = y “^T V s 


( 3 . 4 ) 


where £2 = curl q, and the specific entropy s has been made ziondimensinn- 
al with the gas constant R. Since the upstream value of V contributes terms 
of higher order than those to be retained here, the shock-polar equation 
becomes 


= (U -U,)^ 
d u d 


U U, -1 
u d 


2uy(v +1) . (U^U^ - 1) 


( 3 . 5 ) 


where the subscripts u and d here denote, respectively, values immediately 
upstream and downstream of the shock wave. Since the jump in the velocity 
vector across a shock wave is in a direction normal to the shock, the shock-, 
wave slope is 


dX V, 

B_ d 

dY " U -U. 

u d 

where the shock-wave location is denoted by X = X^(Y). 
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If the nondimenaional friction velocity is small in comparison with 
the nondimenslonal shock-*vave strengtli t , the sonic line in the undisturbed 
boundary layer is very close to the wall, as can be seen from Eqn. (2.14) 
and Fig. 1. That is, if ■* 0 and u^/e -♦ 0, then also 6^ /6 -* 0. A com- 
plete description of the local pressure changes would require both an "outer'' 
solution, obtained by taking a limit of the equations witli coordinates Y/6 and 
X/A held fixed, and an "iimer" solution, obtained with and X/ A ^ 

fixed, for suitable choices of A and A^ . The shockwave can extend nearly 
to the wall, as shown in Fig. 2, and so the upstream influence described 
by the inner solution is very small; it is shown later that A^ - 
where 6^/6 = 0(oxp(-x c/u^)) from Eqn, (2.14). For the outer solution, 
therefore, can be taken equal to the undisturbed velocity (2.10) or (2.4). 
The inner solution describes portiirbations about the undisturbed boundary- 
layer flow, while the outer solution describes perturbations about a different 
boundary -layer flow, downstream of tlic shock wave; the two solutions 
should match in a proper asymptotic sense. 

For Y = 0(6) the length scale A m the downstream direction is found 

from Eqn. (3.1) and the vorticity equation to be A = 0(b 6), where 
2 2 

b (« ) = 1 - M and M is the Mach number in the external flow behind a 
o o o 

normal shock wave. Coordinates x and y arc defined by 


X 




(3.7) 


, whore 

b^ - (y + l)^/^c {1 - ^ (2y + De + . . . } 


(3.8) 
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Since the shock wave is nearly normal, the shock-polar equation (3. 5) gives 
= (1 + c ) ^ This result suggests that throughout the flow down- 

stream of the shock wave U should be represented as a constant value 
1 - € + , . . plus small perturbations of order u_^. It is convenient to sepa- 
rate the rotational part, which can be calculated from Crocco* s theorem 
(3.4), and the irrotational part, which is to be found from the solution of 
Eqn. (3.1) satisfying the appropriate boundary conditions. In the limit as 
— 0 with X and y held fixed, the velocity components are then expressed 
in the form 

"1 ^ 7*^ 2 i 

U = (1 + e) + u^uj '(x, y;e ) + u^u^ '(x, y;e ) + . . . 

2 

) + •• • (3.9) 

V/bo(€ ) = y;e ) + u^<|> 2 y(x, y;e ) + . . . (3. 10) 

where the functions of € shown will be expanded below for e -• 0, 

The entropy s is nearly constant along a streamline, and the equation 
of state gives P = p T, since changes in p ' T* are of higher order than 
terms to be retained here. It follows that along a streamline P T~^ 

-y / (y 2 

"* to the order required here. Substitution of T = a from Eqn. 

(3. 2) then gives the pressure as 


P/P = 1 - V (U-U ) - y ^(U -J)(U-U ) + ,. . . 
e ' u a u 


(3.11) 


Also, Crocco' s theorem gives J2 ~ y ds/dijj, where i}j ^ = p U, - pV, 

A A 

and n = - Uv' so J2/P ~ fl yp^ along a streamline. Substitution in 


u e 


the expression for fi allows calculation of terms in the rotational part of U; 
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(3.12) 


= (1 + 2\e + ...)uQj^(y) 

_ oo 

“2 = f ^ ■ * ‘ (3.13) 

y 


Substitution of the representations (3.9) and (3. 10) into Eqn. (3. 1) leads to 
differential equations for <}>^ and ({j^: 


<}>, +<|>, =0 ( 3 . 14 ) 

^Ixx lyy 

«{>, +<t>y = -4'"’ + (l+€ ){M^/b )^(2 +(v -l)M^)(u^/^ + (|», )<{,, 

2 XX Zyy Zx o' o o 1 lx Ixx 

Expansion of the shock-wave slope (3.6) gives, after integration, the shock- 
wave location x = x (y;u , e ) as 

6 T 

“t 1 

3=s = u^x^^(y;e ) + . . . = ~ (1 + ~ € + • . . ){‘}>^(0, y;e ) - 4>^(0, l;c )} 


+ ... (3.16) 

where the origin of coordinates has been chosen so that x = 0 at y = 1. 

s 

Thus the shock wave is located at x = 0 in a first approximation, as implied 
in Fig. 2, and the flow properties are to be studied in the quarter-plane 
X > 0, y > 0. Boundary conditions at x = 0 are found from the shock-polar 
equation (3. 5), expanded in Taylor series about x = 0: 

^^^(0, y) = -2{ 1 + (y - l)e + ...}u^^(y) (3.17) 

■I’zJO’y) = -f«)'l>ly(0,y) 


3Lsi 

2 



(3.18) 


It is also required that <{» 
approach zero as x^ + y^ 


ly 


(x, 0) 
' 00. 


= «j>2y(X,0) 


= 0 and that all disturbances 
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The limiting form of <|»j^ as € *♦0 was first given in Ref. [8] ; here a 
term proportional to 6 is included. The solution is expressed in terms of a 
distribution of sources along the y^axis: 


4>^(x, y) = + ...} / UQj^(n)ln{x + (y-rj) 

••oo 

(3.19) 

The extended definition ~ gives a potential for - oo < y < oo 

which is symmetric about y - 0 and thus satisfies the boundary condition 
2 2 

there. As x + y 0, the contribution to the complex velocity is 

u^((|) j^^-i<|> = - 2u^{l +(V -1)€ + . . . }fc ^(In a ~ 2H ) + . . . (3.20) 

where z = x + iy. The pressure at the wall found from evaluation of Eqn. 
(3. 11) as y -* 0 is 


P -P, 

w £ 
P 


= 2y u^{l + (2y -l)c + ...}“ J 


2 2 
X +n 


f ... 


where PVP = 1 + y {2e + (2y -1)€ + . ..} is the pressure ratio across a 

X 6 

normal shock wave when the upstream speed is = 1 + «. At larger dis- 

* 2 . 2 _ 
tances, as x + y oo. 


1 


where m^ is defined by Eqn. (2.12); substitution of the approximate analyti- 
cal form (2. 13) gives = 8(1 + II)/k. That is, the integrated effect is 
that of a concentrated source having nondimensional volume strength per 
unit length equal to {1 + . . . }mj^u^6 . One -fourth of this fluid appears to be 
added to the external flow in the quadrant x > 0, y > 0. Since 


16 



d(p U) ~ (1 - M^)dU along a streamline downstream of the shock wave, and 
2 

1 - M (y +1)c , the local increase in the boundary-layer displacement 
o 

thickness is ^Y m.u 6 + .. ., as can also be found by direct calcula- 
4 XT 

tion. An equivalent observation was made for € = O(u^) by Melnik and 
Grossman [5, 6] . ' Perturbations in turbulent stresses contribute only a 
higher-order chcuige locally; the present result does not include the fur- 
ther displacement effect which occurs on a larger length scale as a new 
equilibrium velocity profile is approached. Finally, the shock-wave shape 
found from Eqn. (3. l6) is, for y 0, 

u . 

Xg(y) - = jj" {1 + (V - y + ... } (3. 23) 


and, for y "*■ oo. 


u m 

x^(y)-*^(0) = U + (v +...){ — Iny 

4 °° 

+ “ / +...}. 
^0 


(3. 24) 


For y = Y/6 -*■ oo, the shock-wave displacement continues to increase, and 
should be matched with a suitable perturbed external-flow solution evaluated 
as Y 0. 

The solution for can be found in two parts. A particular solution 
of the differential equation (3. 15) can be made to satisfy homogeneous 
boundary conditions - 0 if sources are distributed 

over the entire x, y plane with the source strength chosen to be an even 
function of both x and y. The boundary condition (3. Id) at x = 0 is then 
satisfied by a distribution of sources along the y-axis, with strength taken 
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to be an even function of y so that 4 * 2 y I'ennains zero at y = 0, as in the 

solution for OX special interest is the total source strength found by 

2 2 

carrying out tlie integrations as x + y ■* oo, with the help of integrations 

by parts and Eqns. (3.14) and (3. 17). The pressure, correct to order 

2 , 2 2 ,- 1/2 2 2 

u_^(x -f y ) as X + y oo, and the second-order source strength 
are found to be 


P - P 


p — ^ = - V (u^[l + (2 y - 1) e + . . . ] + u^(l 4- . . .)m 2 + ■ . • } 


1 X 

2tt 2,2 
X +y 


(3.25) 


m, = 2(5y +9)/ a?,(y)dy - 2 (y +1) / <}>^ (x, 0)dy 


(3. 26) 


For a constant value of y such that y » li P initially decreases as x in- 
creases from zero, reaches a minimum at x = y, and then increases again. 
However, there is a small error at the shock wave x = x (y) = 0(c In y), 

S T 

2 Z 

because the largest term in Eqn. (3.25) is 0(e In y/y ), whereas the 

correct first approximation is found from the shock-polar equation as 
1-1222 

- “V c u^mj^/(2Try) . If it is desired, the accuracy of Eqn. (3.25) can 

2 2 2 2 

be improved near x = x by addition of a term - y u y (0# y)/(x + y ) 

2 

with y approximated by its leading terms 0(ln y) and 0(1) as 

y -► oa; away from x = x the added term is smaller than the second-order 

term originally shown, 

2 2 

As X + y “*•0, the perturbation velocity becomes large, and it is 


again clear that an inner solution is required. For the choice of origin 
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shown in Eqn. (3. 16) and in Fig. Z, x (0) ^ 0 and so the singularity in Eqn. 

(3. 20) is displaced from its correct location through a distance 

- X (0) = 0(u /e ). The domain of validity near x = 0 can be extended slightly 
a T 

by addition of a term - 2u >c ^(1 + . . .)ln(l - x (0)/z) in Eqn, (3.20) for 

accomplished formally by taking a limit as x 0 
with € x/u^ held fixed and then constructing a composite solution. The cor- 

-1 z 

rection is local, and introduces only a smaller change of order e u_^/|z| 

when [z/x (0) | » 1. The modification is, however, necessary for matching 
s 

with the inner solution. A discussion of the inner solution given in Ref. [8] 
is briefly reviewed here, in a slightly n, ified form. For Y - 0(6^), the 
undisturbed velocity is = 1 + O(u^), and the differential equations show 
that changes in U along a streamline are also O(u^) in a distance 


l/2 

AX = 0{u_^^ Inner variables x=** and y* and disturbance velocities 


0 =** and v* are defined by 


X* = 


o s 


(k X - b^6xj0)] 

1/2 


K T 




(U - 1), 


r* - 






V* = 


(3.27) 


(3.28) 


CV 

T 

1/2 1/2 

where factors (k ) ' have been included for convenience. Equation 
(3. 1) and the vorticity equation are then approximated by the transonic small- 
disturbance equations with prescribed vorticity: 


u*8u*/8x* -8v*/3y* + ... = 0 


(3.29) 
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(3. 30) 


o V 


/9 X 


* , « 
- 9 u /d y 


, * 

1/y + 


Xlie shock-wave relations (3,5) and (3.6) become 
2 

I*!' ^2* 

'^d = I<“a -“d> <“d +“u> +••• 


dx 


_s^ 

* 


dy 


U -U , 

u d 


(3.31) 


(3.32) 


where the subscripts a and d again refer to quantities immediately upstream 

/ 

and downstream of the shock wave and the shock-^ave location is given by 


X* = X (y*). As x’*' -*■ - 00 , u^f* approaches the undisturbed form 
s 

u* In y'5' ; the boundary condition at the wall is v* (x^*, 0) = 0; and as 
X* -*■ CO, y*** -♦ oo the solution should agree with the outer solution evalu- 
ated for X - X (0) 0, y -► 0. 

s 

Although complete solutions for u’^' and v^ can only be obtained numer- 
ically, the asymptotic behavior is found relatively easily upstream as 
x=5* - oo and downstream as x=^* -*• oo, -* oo. As x* *♦ - oo, the solution 

has the form 


■ 1 ^ 

~ In y*f= + e^* £(y’*^), v* ~k f (y ) (3.33) 

2 

where f" - (In y*)k f = 0 subject to the conditions that f' (0) = 0 and that 

incoming disturbances be absent as y* -♦ oo; the latter implies 

1/2 

f‘ ~ - (In y*) ' k f as y* -► oo. Numerical integration gives k = 0. 59. 
Downstream a suitable class of intermediate limits should be studied. As 
-* oo, a shock wave is present and must arproach the nearly normal 
shock wave described by the outer solution. Thus, for oo, since 
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« * if ^ ■ff 

In y , Eqns. (3. 31) and (3. 33) give u (0, y ) ~ - In y . If an inter- 
mediate variable ~ V /t}(a^»« ) is introduced, with 1« rj(u^,e) « 6/6 , 
then In y In rj + In y^, where the first term is large and constant where- 
as the second term is 0(1) and variable. In each of the differential equa- 
tions (3.29) and (3.30) the two largest terms remain of the same order if 
X = 0(n/ln rj ) and v = 0(/in n ); then x and y (In y ) ' are of the 


same order. In the limit as x -* oo and y -*• or> with 


X /(y /In y ) held 
fixed, In y In x and so also x /(y /In y ) ~ x /(y /In x ). For the 
derivation of higher-order terms, not to be shown here, it is convenient to 

# if 

make this replacement. In this limit, then, the largest terms in u + In y 

"S' -I /2 * if . if J 5 

and (In X ) v can be written as functions of x /(y /In x ), The solu- 
tions are easily obtained and the results for U and V finally can be rewritten 


as 


1 /? -1 ♦ 1 /7 -1 !}! -1 »!«* 

U - 1 + (/£ u^ln y - (it ) u^ln {x (In x ) + y } (3. 34) 

V - {(y+l)(K x*}^'^^2(k T^/^)"^u tan"^{y* (In x*)^/^/x* )(3. 35) 

Factors (1 - ~ {(■y+l)(/c In x where M is the local 

© * 

Mach number, appear in the locations expected for solutions of the Prandtl- 
Glauert equation. The flow is represented by superposition of a known rota- 
tional flow and an initially unknown irrotational flow, described in terms of 

z -1/2 

perturbation velocities U - 1 and (1 - M ) V which are linear in 

(kT^^^) and are functions of variables (1-M^) - b^6x^(0)]/6 

2 2 

and Y/6 jjj . For a limit such that [x - x^(0)] + y -►0 sufficiently slowly. 
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2 l/2 

with (1 - M ) ^ y/[x-x (0)] held fixed, the largest terms obtained if Eqns. 

(3. 34) and (3. 35) are rewritten in the outer variables x and y are identical 
to the largest terms found from Eqns. (3. 9) and {3. 10), with the help of 
Eqn. (3. 12) and a modified Eqn. (3. 20) in which z is replaced by x x (0) 

+ iy. Introduction of the inner solution thus removes the logarithmic singu- 
larity which appears in the outer solutions for the velocity and the pressure 
as X, y 0. 

In the derivation of these results, terms 9 (p U' U' )/9X, 9 (p U' V )/d Y, 
etc., were omitted from the momentum equation, and therefore a correspond- 
ing set of terms was omitted in Eqn. (3. 1) and in the calculation of the 
changes in vorticity. Expressions for these quantities, and therefore also 
order-of-magnitude error estimates for the solutions given above, can be 
obtained from the Reynolds-stress transport equations [25] . The equation 
for 9 (pU* U* )/9X contains, in particular, terms proportional to 

p U' U' 9 U/9 X. Ahead of the shock wave p U' U' is expected to be of the 
2 

same order as p U' V’ , of order u^. Relative changes at the shock wave 

have been estimated [26] to be proportional to the shock wave strength and 

2 

are therefore small. Thus, p U' U* = O(u^) downstream of the shock also; 

since 9U/9X - for X = 0(€^^^u^), the product is 0(u^e 

Other terms involving velocity correlations are likewise at most 

Neglected terms in the expansion of Eqn. (3. 1) and the vorticity 

T 

equation are also of this order, and can easily be shovm to be small in 
comparison with any of the terms retained. Similarly, the derivative of 
the entropy along a mean streamline contains terms proportional to 
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p U' V 9U/9 Y, etc., and therefore is small enough to be neglected in the 
derivations above. At the shock wave the entropy jump for Y = 0(6) con- 
tains a constant term of order and functions of y which are of order 
2 2 

e * ' * * then be shown that these changes are also suf- 
ficiently small that £2/P and PT remain constant along a mean 

streamline to the order considered here. Finally, the chcuiges in Reynolds 
stresses become important in a sublayer where the perturbation in is 
no longer negligible in comparison with the perturbation in p UU^. For 

X = 0(e^^^u ), since r = O(u^) and U^- = 0(e the sublayer is de- 

2 1 /2 

fined by Y = 0(u^c * ). As will be shown in detail in Part II, the relative 
change in r is 0(e), and the new term in U which contributes to a displace- 
ment effect is 0(e u_^). From the continuity equation it follows that the 

2 2 

corresponding term in V is 0(e u^). Thus, as y ■* 0, the largest term 

in the outer solution for V which satisfies a nonzero boundary condition 
2 2 

is 0(e u^), smaller than any of the terms retained above. All of the 

neglected terms arising from these effects are smaller than the terms 
retained by at least a factor of order e. 
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4. Geometric Effects and Comparison with Experiment 

The theory of the preceding section leads to a limiting form for the 
pressure distribution as u^ -*■ 0 and 0# for unseparated flow. In the 

flow past an airfoil at supercritical speed, with a shock wave terminating 
a region of supersonic flow, the additional effect of surface curvature can . 
also be important in changing the pressure distribution and delaying sepa- 
ration, as discussed below and in Part II, The boundary layer might remain 

attached for M up to about 1. 25, depending on the profile shajpe; Re may 
7 8 

be about 5 x 10 or perhaps as high as 10 ; and the flow ahead of the shock 
wave, experiences a favorable pressure gradient, with magnitude which de- 
pends on the airfoil shape, so that the profile parameter II is 

smaller than 0. 5 (e, g,. Ref. [30] ). For a combination of parameters which 

is favorable with regard to requirements of the present theory, with M = 

c 

g 

1, 26, Re = 10 , and II = 0, the relative position of the sonic line is given 
by 6.^ /6 =0. 10. This value would increase as or Re decreases or as 
H increases, as seen from Eqn. (2.14). Experimental results, however, 
are not yet available with detailed local pressure measurements for values 
of the parameters which correspond to such airfoil flows and which meet 
the requirements of the theory. For all available data, either the flow is 
separated or the values of the parameters are such that the sonic line is 
not close to the wall. Nonetheless, a comparison with data from Refs, [l] 
and [2] has been carried out, and the agreement seems favo rable provided 
that corrections for geometric effects are included. 
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A wall having convex longitudinal curvature is described locally by 
1 2 

Y ~ ” KX , where K « 1 if the radius of curvature is large in comparison 

with the reference length used in the definitions of X and Y. A local solu- 
tion for the inviscid external flow near the foot of a normal whock wave shows 
a discontinuity in streamline curvature [27, 28] . Ahead of the shock wave 
Py > 0 to provide the required acceleration toward the wall; if the flow is 
irootational, it follows that Uy < 0, The shock-wave relations give Uy > 0 
and Py < 0 downstream; therefore also > 0, whereas the tangency con- 
dition at the wall requires < 0 as Y-»- 0, The term in the complex 
velocity which satisfies the required conditions as X, Y -► 0 has the derivative 

’^o^X^ " ^^X^ ~ In Z + iK + C(K) (4. 1) 

1 -1 2 
for 0 < arg Z < “v, where Z = b^ X + iY and, as before, b^(®) ~ (\+l)^» 

The largest omitted term is of order K and is real; the value depends on 
the flow description for Z = 0(1), and is known for symmetric two-dimen- 
sional or cixisymmetric nozzle flows [29] • 

Terms and of order K€ u_^ and Ku^, respectively 

when Z = O(u^), are now added to the expansions of U and V given by £qns. 
(3.9) and (3. 10). The rotational part of U is unchanged, and reformulation 
of the boundary -value problem for the perturbation potential shows that 
is unchanged, whereas now depends on K, through nonlinear terms in 

the potential equation; that is, <|»2 " new terms in 

2 

contribute a change in U which is O(Ku^), smaller than terms retained pre- 
viously provided that K = p(l). Thus, to the order considered here, for 
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u^/e -*■ 0, a cur va tare correction is simply added to the earlier results. 

The new term in the pressure, written in terms of x and y, is 

{ 2x In 6 + X ln(x^ + y^) + T y - 2y tan ^ ^ + Ax} (4. 2) 

where the constant A is determined only if a solution is known for the external 
flow at larger distances. 

An early careful and comprehensive expe-^imental study was carried 
out by Ackeret, Feldmann, and Rott [l] . In Figs. 3 and 4, predicted pres- 
sures are compared with their experimental results for = 1.32, corre- 

5 

sponding to € = 0. 247, and Re = 9. 6 x 10 , based on distance to the shock 
wave. Eqns, (2.9), (2. 11) and (2, 14) are used for approximate evaluation 
of other parameters. One more experimental value is needed; is chosen 
since it is easily read from the measured velocity profile and since oaixLy 
In 6 ^ enters the equations, so that an error has small effects on other quan<^ 
titles. For 6^ = 0.0055, the calculations give u^ = 0.051, 6 = 0.021, and 
n = 0.28. This value of II seems plausible (e.g. , Ref. [30] ) because of 
the observed small favorable pressure gradient ahead of the shockwave. An 
adverse gradient of about the same magnitude is evident downstream, and is 
estimated here by P ^9P/3X a 0.12, where P is the upstream stagnation 

t c 

pressure. A corresponding term is added to Eqn. {4. 2) and the term propor- ‘ 
tional to K6Ax is neglected. The local curvature of the plate can be inferred 
from measured pressures immediately behind the shock wave. It is estima- 
ted that P^^9P/9 Y a 0.15; since Py ~ it follows that K a 0,2, With 

the kind assistance of Prof. Z. Plaskowsfci of the InstitQt fiir Aerodynamik, 
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ETH Zurich, the author was able to measure ordinates of the plate actu- 
ally used in the experiments: values in an appropriate neighborhood confirm 
tVie estimate K a 0.2. The origin x = 0 is chosen at the estimated position 
of the shock wave at the edge of the boundary layer, found using measured 
pressures outside tlie boundary layer together with Eqn. (3.24). 

The comparison in Figs. 3 and 4 shows that the curvature effect is 
comparable in importance with the boundary-layer displacement effect; 
addition of the curvature term leads to a more pronounced "shoulder" in 
the predicted wail pressure distribution. Tlie longitudinal pressure grad- 
ient due to tuiinel divergence is also seen to be important. At the plate for 
typical values of -x, say 4 < x < 14, the prediction gives about 75 percent 
of the pressure drop below the value for a one -dimensional flow; outside 
the boundary layer, at Y/6 ~ 3.6, the agreement is somewhat better. It 
is found that the velocity in Eqn. (3.21) is closely approximated by (const, )/x 
for X > 2, so that Eqn. (3.25) for the pressure is adequate here, with the 
correction (4.2). Modest changes in the assumed values of the parameters 
do not have a major effect on the comparison; for example, at a giveii X, 
m^^/x does not depend strongly on n because 6 increases if II decreases. 

The upstream exponential decay predicted by Eqn. (3. 33) is also shown in Fig. 3, 
in the form A p/p ~ u exp{k,(x -x )}, with x taken equal to -14 for 

approximate agreement with experiment. A major difficulty with this com- 
parison is that the upstream sonic line lies at about y = 6,j,/6 = 0.26, and 
the shock wave ends at a still larger distance from the wall, so tliat the 

* <c 

inner region for x = 0(1), y = 0(1) is not negligibly small. At a higher 
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Reynolds number and therefore a lower the shock wave would extend 
closer to tlie wall, and the size of the region in Fig. 3 where no prediction 
is given would be smaller, A second serious difficulty arises because the 
flow probably was separated. The authors of Ref, [l] stated that reversed 
flow would not be ascertained at any point; however, the velocity profiles 
shown seem inconclusive, since measurements were not possible very 
close to the wall. Calculations based on the theory of Part II of the present 
paper, for the parameter values given in Ref. [l] , indicate that the flow 
was in fact separated, with a very thin separation bubble having length 
equal to a few boundary-layer thicknesses. The effect of such a bubble 
would give a more gradual pressure rise in the region of greatest disa- 
greement in Fig. 3. Finally, a slight unsteadiness in the shock wave posi- 
tion would also contribute to a decrease in the measured pressure gradient. 

A correction for flow in a circular pipe can be derived in terms of 

4 " 4 “ 

cylindrical coordinates x and r defined by 


+ 

X = 


b (€)R 
o 




1 Y _ , 6. 

" = ^ 'R - ^“Ry 


(4.3) 


where R is the ratio of the local pipe radius to the reference length, and Y 

4 " 

is measured inward from the wall, so that r = 0 at the axis. Solutions 

+ + 

are to be found for c -*■ 0, 0 with x , r fixed. The wall shape is 

+ 2 + 
given by r = 1 + c f(X/R) with f = 0 at X/R = 0. Velocity components U , 

4 * 

V in the x , r directions can be written with the local curvature and 


boundary-layer effects shown separately: 
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= U^‘^^X/R, r\i) + r^) + ^ ^ ^ t r'*') 


+ . . . 


(4.4) 


= V^°\x/R, r^c) +Kb^(€)V^^^x’*', r'*’) +a^|b^(€)V*^^x^r'*') +... 

(4. 5) 

2 

where now K = e f’(0) is the wall curvature at the foot of the shockwave, 
made nondin:.ensional with the reciprocal of the pipe radius. The terms 
and are the terms which would be present if the effects of the 

shock wave were ignored [29] . Terms proportional to K contain the local 
curvature effect, and terms proportional to u^6/R contain the local boun- 
dary-layer displacement effect. The latter is described in terms of a ring 
source of radius r =1 located at x =0 and having volume strength per 
unit length equal to mj^u^6 + . . , ; numerical solutions for ~ 0(1) given 

by Melnik and Grossman [6] also include this effect. For x -* oo, the fluid 
added at the source gives an increase of “m^u^6/R in U^, shown explicitly 
in Eqn. (4. 4). 

l/2 

The local solutions for X/R = 0(c ‘ ) are found in terms of a stream 
function defined by where i=l, 2. 

The largest terms in Eqn. (3. 1), combined with the irrotationality condition, 
tion, lead finally to 


CO -X X 

I (0/ + +x y' (i) n +T i\ \ 
iL (x,r) = Z-/a e rJ{\r) 

, n in 

n=l 


(4.6) 
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where ^ • • • » that the wall boundary condition 

1) = 0 is satisfied; also -*• 0 as oo, and boundary values 

are to be specified at = 0, To the order required, the shock-polar equa- 

(1) ■l‘2 

tion reduces to the Prandtl relation, and so tjj = r (1 - r )/4 at 

= 0. The condition that the ring source gives no term of order u 6/R 

Z ^ 

in at x^ = 0 implies ijj = - m^r /4 at x^ = 0, Comparison with the 

( 2 ) 

wall boundary condition shows that i}i is discontinuous at the foot of the 
shock wave x = 0, r =1; the value obtained as x -► 0, r -*■ 1 depends on 
the direction of approach. The coefficients a^^^ can be found from the solu- 
tions of Messiter and Adamson [29] or by direct calculation, and the coef- 
( 2 ) 

ficients are found directly: 


( 1 ) _ 


n 


\ J ) 
no n 


( 2 ) 


m. 


n 2\ J (X ) 
non 


(4.7) 


for n=l, 2, 3, . . . , 

For calculation of the pressure distribution and the shock-wave shape, 

(2) 

it is convenient to introduce the corresponding velocity potential <)> , which 

satisfies 9 and 9 and which has a loga- 

4 " 4 " 

rithmic singularity at x = 0, r =1. With the help of the asymptotic form 


for J (\ r ), one can show the singular part explicitly: 
o h 
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- - 


m, oo J (\ r ) 

f S (f 

n=l n J (X ) 


-\ X 

n 


jn 


IzH 


n 


7r(n + ^)(r'^)^/^ 


-7T (n + ~)x 

e 


+ 

cos[7T(n + ^)r'‘' - |-1 } 


- ITT? £ w - 4? } H. 8) 

2 ;r(rV/^ 

where In = - (^ /4){x^ + i(l - r^)}, and indicates that the real part 

+ + 

is to be taken. As x -• 0 and r -►1, the largest term in the complex 

(Z) iZ) 

velocity U' - iV is due to a two-dimensional source of strength m^, 

in agreement W'ith Eqn. (3.22). The change in the boundary -layer displace- 

(2) (2) 

ment effect is then found by subtracting the source term from U - iV 

■f* 

and adding the constant term which remains as x ■* oo. If the numerically 
small contribution of the infinite series is omitted, the corresponding 
correction to the wall pressure is 


AP 
w 

P 


e 


- y 



m 

1 


2 



“Stt x^/4 


1 


+ 

-7T X 
e 



(4.9) 


+ • 1 / 
As X = X/(b^R) "*■ 0, AP^/P^ approaches a constant value - ~-y m^a^6/R, 

which implies an additional second-order correction to the boundary-layer 

solution found in Section 3 for X = 0(b 6). The shock-wave shape is found 

o 

directly from the potential; in particular, as r^ 1 the displacement of 
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the shock, from its intersection with the axis is found by adding the pertur- 
bation potential from Section 3 to that found here, and subtracting the com- 
mon term proportional to In Y. The result is 

AXg = - {47re)“^(Y+l)^/^6^/^mj^u^6 ln(R/6) + 0{u^6c“^/^) (4.10) 

In Fig. 5 a comparison is made with pressures measured in a circu- 
lar pipe by Gadd [2] , for M = 1. 12 and c = 0.097. The length of an equiv- 
alent flat-plate boundary layer is not a given quantity; instead, Gadd' s 
estimated value for boundary-layer thickness is used here, along with the 
estimate 6^/6 =0.45 found from the measured velocity profile. The sonic 
line is therefore still further from the wall than in the Ackeret experiment. 
Other approximate values are calculated as u_^ = 0.04, 6 = 0.02, and 
n =0.1; the Reynolds number corresponding to these values is Re = 

6 X 10^. A pressure gradient due to small divergence of the test section 
is estimated downstream by P^. 8 P/8X = 0. 06. The ei^ect of finite pipe 
radius is seen to be about as large as the boundary-layer displacement 
effect. The upstream exponential decay is also shown, with x^ again taken 
equal to -14. Again there is a relatively large region where no prediction 
is made and where numerical solution of the transonic small-disturbance 
equations is required. Such a solution was obtained by Melnik and Gross - 
man [4] for this case and is also shown in Fig. 5. For large x the analy- 
tical and numerical solutions differ by an amount about equal to the correc- 
tion for the change in pipe cross-section area. 
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5, Concluding Remarks 

The interaction of a turbulent boundary layer with a weak normal shock 
wave has been described here and elsewhere [3, 4, 8] in terms of a rational 
approximation based on systematic asymptotic expansion procedures. The 
interaction is characterized by two small parameters, a nondimensional 
friction velocity u^ and a nondimensional shock-wave strength'^, and limit- 
ing forms of the local solutions can be studied as u_^ -► 0 and e -*• 0. For the 
case u_^/« oo [3] , analytical solutions indicated that separation does not 
occur; solutions for u^/e held fixed [4] , with the first approximation de- 
scribed by the transonic small-disturbance equations, gave the same result. 
If, finally, it appeared that analytical solutions would be possible 

and that perhaps the onset of separation could be discussed. Solutions for 
the pressure have been obtained here and will be used in Part II for the cal- 
culation of wall shear stress and a discussion of incipient separation. 

The largest terms in the pressure, of order u^, are derived quite 
easily, and a number of higher-order effects have been added. Corrections 

of order e u_^ give, ‘e. g. , a 35% change if = 1. 25. A partial solution for 
2 

terms of order u^ shows that these terms likewise are significant, typically 

6 7 

giving changes of 25% to 50% for Re = 10 or 10 , Corrections of order 
Ko.^ ajnd u_^6 /R, obtained in \nalytical form for a wall with longitudinal cur- 
vature and for a circular pipe respectively, are found to be numerically 
important for the tests of Refs, [l] and [2] . 
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In the solutions for these higher^rder terms, the dependence on the 
parameters is of coarse shown explicitly, and the relative Importance of 
different effects is therefore apparent. It is not, however, possible to ob- 
tain analytical solutions in the asymptotically small inner region which 
accounts for the upstream influence. For values of the parameters cor- 
responding to actual transonic flight conditions, it is possible for this re- 
gion to be relatively small. Experimental results, however, are not avail- 
able in this parameter range; for existing data, either the flow is separated 
or the sonic line is not close to the wall. Nonetheless, some comparisons 
with such data were attempted, and the agreement seems fairly good down- 
stream from the inner region. The predicted pressures remain somewhat 
higher than the experimental values, and the correction terms calculated 
thus far are large enough to suggest that additional higher-order terms would 
be likely to give still further improvement. 

An essential feature of the asymptotic flow description in terms of u_^ 
and € is the two-layer structure of the undisturbed profile, expressed by the 
law of the wake and the law of the wall. It is this property which permits 
the calculation of interaction pressures without knowledge of changes in 
shear stresses close to the wall. In other studies [31 , 32] which were not 
based on use of this profile, derivation of a sublayer solution was necessary 
before the calculation of the pressure could be completed; these studies also 
introduced a linearized formulation for the main part of the boundary layer. 
In the present asymptotic description for u^/f -^ 0 , sublayer effects do not 
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appear even among the second -order terms in the pressure, JLinear equa- 
tions appear naturally as a consequence of the limiting case considered, 
and the procedure for adding higher -order terms is clear. In the formula- 
tion for held fixed [4] , again the flow details near the wall do not 

influence the pressure. The differential equations obtained in the limit 
are, however, nonlinear and numerical solution is required. For u A -► cx> 
[3] , the pressure once more is found without knowledge of changes in the 
wall shear stress, 

A complete asymptotic description for 0< u^/e < oo is therefore now 
available, with numerical solutions obtained as u_^ -*• 0 and c 0 if u is 
fixed and analytical solutions if oo or 0. For accurate 

calculations in parameter ranges of practical interest, some further ex- 
tensions appear to be needed. In the present case, as u /e -*■0, the neces- 

T 

sary condition that the sonic line be close to the wall is met for a rela- 
tively narrow range of the parameters.' The solutions for wall pressure would 
be more useful if a simple curve fit were introduced for the inner region, say 

by means of a straight line tangent to the source solution downstream and to 

* 

the exponential solution upstream. The choice x^ « -14 in the exponential 
term was made for agreement with experiment in Fig. 3; the results shown 
in Fig. 5 suggest that the magnitude is too large and that perhaps a more 

a. 

suitable tentative value would be = -10. The present solutions also suggest 
that terms of higher order than those retained in Ref. [4] are likely to be 
important for =0(1)- this case the curvature correction would no 

longer have a simple form, but would have to be incorporated in the numerical 
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solution through the use of modified boundary conditions. Moreover, it ap- 

2 

pears that certain terms of order u^, and possibly still other higher-order 
corrections, will also be essential for numerical accuracy, in the wall 
shear (Part II) as well as in the pressure. Finally, the local interaction 
influences the potential flow at larger distances; the manner of introducing 
corrections in the external flow deserves further study. 
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Figure Titles 

1. Undisturbed velocity profile. 

i. Asymptotic representation of "outer" flow. 

3. Pressure at wall with longitudinal curvature; = 1.322, Re - 

9. 6 X 10^, K a 0. 2. 

4. Pressure outside boundary layer, at Y/6 * 3.6, for curved wall 

M = 1. 322, Re = 9. 6 X 10 , K == 0. 2. 
e 

Pressure at wall of circular pipe; = 1. 12, Re » 6 x 10^, 

6/R ^ 0.055. 
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Uncorrected (1-Dim.) Flow 
Eqn.(3.25),forPloneWoll 
Eqn.(4.12) Added for Pipe 
Experiment [2] 

Numerical Solution [6] 
Pg/Pf+u^ exp {0.59(x*+14)} 
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